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Abstract 

A Monte Carlo program solving Boltzmann equation for partons via 
cascade method is presented. At presented, only gluon-gluon elastic 
scattering is included. The scattering cross section is regulated by a 
medium generated screening mass. Three different geometric modes 
(3 dimension expansion, 1-d expansion, and scattering inside a box) 
are provided for theoretical study of the parton transport and the 
applicability of the cascade method. Space cell division is available to 
save the number of computer operations. This improves the speed of 
the calculation by a large factor and makes the code best optimized for 
simulation of parton cascade in ultrarelativistic heavy ion collisions. 
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Program Summary 



Title of program: ZPC 1.0.1 

Program obtainable from: CPC Program Library, Queen's University of Belfast, 
N. Ireland and from the author at bzhang@ntl.phys.columbia.edu 

Computer for which the program is designed: SPARC stations; Installation: 
Physics Department, Columbia University, USA 

Computers on which it has been tested: (1) IBM RS/6000; (2) DEC Personal 
Workstation 600 AU; Installation: (1) RHIC Computing Facility, Brookhaven 
National Laboratory, USA; (2) Theory Program, Nuclear Science Division, 
Lawrence Berkeley Laboratory, USA 

Operating systems under which the program has been tested: SunOS 4.1.3, 
AIX 1 

Programming language used: FORTRAN 77 

Memory required to execute with typical data: 1,450 words 

No. of bits in a word: 32 

Number of lines in distributed program: 6330 

Keywords: ultrarelativistic heavy ion collisions, partons, Boltzmann equa- 
tion, screening mass, parton cascade, cascade schemes, collision frame, or- 
dering frame, ordering time. 

Nature of physical problem 

In ultrarelativistic heavy ion collisions, perturbative QCD predicts the pro- 
duction of a large number of minijet gluons. The gluon system produced near 
central rapidity can be considered as an almost longitudinal boost invariant 
and baryon free system. Gluons inside this hot matter will interact among 
themselves before hadronization. The effect of the final parton interaction on 
the global variables should be studied to match the experimental observables 
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to the pQCD predictions. 
Method of solution 

The space time evolution of final state partons can be approximated by a 
Boltzmann equation. To solve the Boltzmann equation, we employ the cas- 
cade method in which two partons scatter when their closest distance in a 
specified frame is less than the interaction length. Gluon scattering cross 
section is regulated by a medium generated screening mass. Space cell di- 
vision is made to speed up the calculation. Options to do the cascade with 
and without space cell division are provided to check the applicability of the 
optimization. 

Typical running time 

The running time depends on the density of partons and the scattering cross 
section. For 3 dimensional expansion, when the initial density in the local 
rest frame is ~ 50/ fm 3 (at proper time 0.1 fm/c, with transverse radius 
5 fm), and the parton parton scattering cross section is 3.5 mb, with space 
cell division, it takes around 6 minutes to finish 1 event 7000 binary col- 
lisions) on SPARC20, and 30 seconds on alpha Personal Workstation. 

Unusual features of the program 
none 
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Long Write-Up 



1 Introduction 

In ultrarelativistic heavy ion collisions at collider energies, hard and semihard 
processes dominant. The high energy density produced in these collisions 
opens the possibility of studying many interesting phenomena of strong in- 
teractions, e.g., deconfmement phase transition and chiral symmetry restora- 
tion. The space time evolution of the hot matter produced in the collisions 
may make the final observables drastically |1J different from those directly 
predicted from pQCD and Glauber geometry. So it's very important to study 
the final state interactions in order to link the experimental observables and 
the pQCD predictions. 

The final state parton interactions can be described by the equations 
of motion for the quark and gluon Wigner operators. Under semiclassical 
approximation and the Abelian dominance approximation, in the weak field 
limit, the equations of motion become a set of equations resemble the classical 
transport equations]^]. This motivates the study of final state interaction by 
solving the Boltzmann equation for quarks and gluons: 

P^d^f a (x, p, t) = J2m 2&i,&2,-,&m / ITS 1 (2n)32E b f b i 0^> P*>i ' 

V V l'U n d3pcj I M I 2 

l^n l^ci,c 2 ,-,c n J W.?=l (27r) 3 2_B Cj \ lvl m-*n\ 



(27r) 4 5 4 (E™=i^-Er=iP Ci ) 

• sr=i w 3 (p - Ph) + n=i ^np - pd] ■ (i) 

In the above equation, f a (x, p, t) is the phase space distribution of parton 
type a, and M m ^ n is the scattering amplitude defined by: 



< Ci, C2, ■ • ■ , c n \S\bi, 62, ■ • ■ , b m >- 



(m n 
Y,Ph-T,Pc 
k=l 1=1 



To study the final state parton interactions, we have developed ZPC, 
which solves the Boltzmann equations for partons by the cascade method. 
In the cascade model 0, two partons scatter when their closest distance is 
less than J£ (er is the parton parton scattering cross section). The scattering 



4 



angle is determined by the differential parton parton scattering cross section 
||]. Because cascade prescription of solving the semiclassical transport equa- 
tion doesn't assume any particular form of the phase space distribution, it 
can be used to study systems far from local thermal equilibrium and it takes 
into account of finite mean free path effect automatically. 

In the inside outside cascade picture, the partons come out as indepen- 
dent partons after the nucleus-nucleus collision. These partons are taken as 
the initial condition for the parton cascade. The parton formation proper 
time usually depends on the transverse momentum of the parton ||. In the 
simplest picture, the partons are formed on a hyperbola of constant proper 
time H determined by the average transverse momentum. 

The space time evolution afterwards depends also on the scattering cross 
section. Currently, a rough Yukawa type of screening is used to regulate the 
parton scattering cross section. Lattice gauge theory |7J and semiclassical 
transport theory || can give us screening mass(es) for parton interactions 
for equilibrium system or system close to local thermal equilibrium. Assume 
that the formula can be used for the far from equilibrium case as in the 
ultrarelativistic heavy ion collisions, we can get an estimate for the medium 
generated screen mass. On other hand, this provides one possible way of 
measuring the effective screening in the medium from the experiment. 

In ultrarelativistic heavy ion collisions, there are indications that multi- 
parton processes other than elastic scattering may dominate . Also because 
inelastic collisions are important for the chemical equilibrium, it is important 
to look at parton production and annihilation processes. Besides inelastic 
collisions, parton statistics and hadronization may all affect the prediction 
in a non-negligible way. These elements will be incorporated in ZPC in the 
future. 

ZPC is developed in parallel with GCP by Yang Pang ||. The major 
difference is that ZPC 1.0.1 focuses on gluon elastic scatterings, while GCP 
is being designed to be a more general, user friendly package. ZPC 1.0.1 has 
the local interaction list optimization. ZPC 1.0.1 is written in FORTRAN, 
while GCP is written in C. 

In this paper, we will discuss the physics background for the parton cas- 
cade simulations, some ZPC technical points, then a detailed description of 
the program followed by the instructions for users. 
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2 Cascade simulations of URHICs 



For Ultrarelativistic Heavy Ion Collisions (URHICs) at RHIC energies and 
beyond, the gluon number is much larger than the sea quark and antiquark 
number, and the gluon gluon elastic cross section is much larger than the 
gluon-quark transition cross section ||. To first approximation, we consider 
only gluon gluon elastic scatterings. The initial conditions can be read from 
an external file or some ideal initial conditions can be generated within the 
code. 

2.1 Formation time 

The formation time of partons gives the time beyond which the radiated 
partons can be considered as independent partons. There are different mod- 
els for the formation time. A simple model motivated by the Uncertainty 
Principle and used very often for ideal theoretical calculations is proposed 
by Bjorken||. In this model, partons are formed on a hyperbola of constant 
proper time. The integrated formation probability can be written as: 

P(r)=6(r-r ). 

tq is inversely proportional to < itlt >, or can be simply taken as 1/ < 
rriT >■ Another more realistic model proposed by Gyulassy, Wang, and 
Pliimer ||has an integrated probability: 

P(r) = 1/(1 + (r / /r) 2 ), 

where ry = cosh y/rriT- The above formula interpolated between the per- 
turbative result Par 2 for small times and the asymptotic result which is 
characterized by a formation time cosh y/rriT- 

In ZPC 1.0.1, the ideal initial conditions are generated with a fixed for- 
mation proper time. Options are available to generate formation time for 
parton configurations generated at t — and propagate the partons to their 
formation time. 
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2.2 scattering cross section 

The gluon-gluon elastic scattering differential cross section is: 

da 7ra| 9 / ut us st \ 
~dt = ^ 2 ~2\ ~ 1? ~ ~¥ ~ tf) ' 

At fixed s, the leading divergent term is: 

da 9na 2 s ( 1 I 



dt 2s 2 \t 2 u 2 

Since two outgoing particles are identical particles, the scattering angle goes 
from to |. If we let the scattering angle go from to ir, then the t, u 
singular terms are redundant. Hence, 

da 9na 2 s 1 



dt 2s 2 t 2 ' 

The singularity is regularized by a medium generated screening mass. We 
end up with a scattering differential cross section of the form: 

da 9na 2 s 1 



dt 2s 2 (t-/i) 2 ' 

This is just the Yukawa type of scattering cross section. Currently, we set 
the total cross section to be s independent: 



9na 2 s 

and the differential cross section changes correspondingly to: 

da 97TQ;| ( ^ | fi 2 \ 1 



1 + 



2 ' 



dt 2s 2 \ s J (t-fi) 

This is only a simplified model for the gluon-gluon scatterings. Following the 
instructions for users, it's easy to make the cross section more realistic. 
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2.3 Screening mass 



When the gluon system is not far from equilibrium, perturbative calcula- 
tion relates the electric screening mass m with the phase space density f(k) 
through: 

m 2 = -^flim lq ^ f d'k^lq- Vj/tf). 
71" J q ■ k 

Assuming the formula can be used for nonequilibrium as well, the phase 
space density f(k) is related to the gluon invariant momentum distribution 
g(y, k T ) = d 2 N G / ' dydk\ through: 

9gv \k\ 

At collider energies, the transverse screening mass is almost the same as the 
longitudinal screening mass. In the ideal case that the rapidity distribution 
is an inverted square well and the transverse momentum distribution is an 
exponential, the screening mass can be estimated through: 

mT * 2F- 

With a ~ 0.4, ~ 300, which is typical for the minijet gluon system 
produced at RHIC, we get the screen mass fi ~ 4.5/m~ 1 . So in the cascade 
calculation, the screen mass is taken as several inverse fermis. 



3 General features of ZPC 1.0.1 

ZPC 1.0.1 provides the user with a choice of several scattering prescriptions, 
the option to read in initial conditions or to generate some ideal initial con- 
ditions, and the space cell division optimization. 



3.1 Scattering prescription 

Because of the geometric interpretation of the scattering cross section, in- 
formation travels across the distance of closest approach when a scattering 
happens in a time period of the fixed time step or the mean free time. This 
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may lead to causality violation. Different prescriptions may give different 



results 10 



To specify the scattering scheme, we need to specify the collision frame, 
the collision space-time point, the ordering frame, the ordering time. The 
collision frame is the frame in which the collision happens. The ordering 
frame is the frame in which the collisions are ordered in time. It is the global 
frame where the observer is in. The ordering time is the time of a collision 
that we use for ordering the collisions. 

In ZPC 1.0.1, the collision frame and space time point choice is set 
through the first digit of the variable iordjsch which can be specified in 
input file zpc.ini. The second digit of iordjsch gives the ordering time. The 
ordering frame can be chosen by setting the variable ibst -flag. 



3.2 Geometric configurations and built-in initial con- 
ditions 

ZPC 1.0.1 provides 3 different geometric configurations to do cascade. We 
can study 3 dimensional expansion, 1 dimensional expansion and particles in 
a box by setting the variable Lconfig. Setting igen_flag gives some built-in 
initial conditions. 

For the 3-d expansion, the built-in initial condition is like the minijet 
gluon system generated in RHIC Au on Au central collisions. 4000 particles 
are uniformly distributed in the pseudo-rapidity range —5 to +5. They are 
generated in local thermal equilibrium with a temperature 500 MeV (h = 
c = 1). At t = they are within a transverse disk of radius 5 fm. The 
particles are formed with a longitudinal formation time r = 0.1 fm. 

The 1 dimensional expansion built-in initial condition is different from the 
3-d one in that the transverse positions of particles at their formation times 
are within a square with x (y) goes from —5 x sizel (2) to +5 x sizel (2). 

The particles in a box built in initial conditions give particles in thermal 
equilibrium with temperature 500 Mev. The particles are confined within a 
box with x (y,z) goes from —5 x sizel (2, 3) to +5 x sizel (2, 3). 



3.3 Optimization for fast cascade code 

The parton cascade is different from the hadronic cascade for relatively low 
energy nuclear physics in that there are much more particles to start with 
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in one event. For example, at RHIC energies, each central Au on Au event 
has around 4000 minijet gluons. To get the next collision partner for a 
particle after the collision, a number of N (N is the total number of particles) 
operations have to be performed. The total number of operations for one 
event is proportional to N 2 . For N ph 4000, it's difficult to get enough 



statistics for detailed study of events, e.g., correlations [11]. To improve the 



speed of the parton cascade, we chose to divide the space into cells. In this 
manner, to update a particle's collision partner, only the particles in the 
same and neighboring cells need to be checked. This makes the total number 
of operations go down to C*N*n (n is the average number of particles in one 
cell). 

Cells expand with time after time size (which is the minimum of the cell 
dimensions sizel, size2, and size3). But to ensure the results are exactly 
the same as the results without cell division, size should be larger than the 
interaction length. For reference, the user may want to produce the results 
without cell division for reference and cross check. For the 3-d expansion 
case, set v 1 = v2 = v3 = sizel = sizel = size3 = 0; for 1-d 
expansion, set i-config = 3; for particles in a box set i-config = 5. 



3.4 Optional OSCAR standard output 

OSCAR (Open Standard Codes and Subroutines) QT2J standard has been 
proposed by Yang pang and passed by a group of physicists (Open Standard 
working Group) in the field of cascade simulation of high energy nucleus- 
nucleus collisions to provide objective scientific criteria for cascade simu- 
lations. Currently, the standard output (OSCAR1997A) has a file-header 
and an event-header. The file-header contains general code and simulation 
information and the event header stores particle number, particle ID, 4 mo- 
mentum, mass, position and time. We can set i-oscar = 1 from file zpc.ini 
to write out the standard output into the file zpc.oscar. This provides a way 
of reconstructing event history. 



4 Program description 
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4.1 Main subroutines 

To clearly analysis the program structure, several technical words must be 
defined in the ZPC context. An EVENT means the time evolution of a par- 
ticular initial phase space distribution of particles. For one particular event, 
the user can test the effect of different random number sequences by running 
the same initial conditions with different random number sequences. The dif- 
ferent space time evolutions with one same initial phase space distribution, 
different random number sequences are called different RUNs. An OPERA- 
TION is one of the following: the formation of a particle, the collision of a 
particle with cell boundary, the collision of two particles or a combination of 
these simple operations. PARTICLE INDEX is the label of particle in the 
common block /prec2/ and PARTICLE NUMBER is the particle label in the 
common block /precl/. The following main subroutines are all in the main 
program zpc and control the flow of the program. 

• SUBROUTINE INLZPC 

Purpose: to initialize the program. It reads in the parameters from the 
initialization file zpc.ini and initializes the main common blocks. 

• SUBROUTINE INLEVT 

Purpose: to initialize one particular event. It reads in initial conditions 
for an event and, or generates initial conditions depending on the user's 
requirement. Then it boosts the initial conditions if the user requires. 

• SUBROUTINE INLRUN 

Purpose: to reset the initial conditions to the same as last run. 

• SUBROUTINE ZPC_RUN 

Purpose: to select the next operation according to ordering time, per- 
form the operation and update the interaction list. 

• SUBROUTINE ZPC_ANA1 

Purpose: to do analysis after each operation. This enables us to look 
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at the time development of physical observables. 



SUBROUTINE ZPC_ANA2 
Purpose: to analyze event final states. 

SUBROUTINE ZPC_OUT 

Purpose: to average results over events and write the results into out- 
put files. 



4.2 Random number generator(s) 

Default random number generator in ZPC-1.0.1 is the linear congruential 



method random number generator rani fl3fl . Another random number gen 



erator ran3 based on the subtractive method is also included. 
4.3 Main common blocks 

There are seven common blocks in ZPC-0.2.1 that contain parameters for 
the cascade simulation. They are named para*, where * is a number going 
from 1 to 7. 

• COMMON /PARA1 / MUL 

MUL: the multiplicity of one particular event. 

• COMMON /PARA2 / XMP, XMU, ALPHA, RTS_CUT2, CUTOFF2 
XMP: the mass of the particles. In ZPC 1.0.1, there is only one type 
of particle, the gluon. 

XMU: the screening mass. 

ALPHA: the strong interaction coupling constant. 
RTS_CUT2 : lower square root of s cut for parton collisions. 
CUTOFF2: the total parton scattering cross section divided by 7r. 

• COMMON /PARA3/ N_STA_EVT, N_EVENT, N_RUN, LEVENT, 
LRUN 

N_STA_EVT: the starting event number in the input file. 
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N_EVENT: the total number of events. 
N_RUN: the total number of runs per event. 

I _EVENT: the current event number that goes from 1 to N_EVENT. 
LRUN: the current run number. 

COMMON /PARA4/ IFT_FLAG, IRE_FLAG, IGEN_FLAG, IBST_FLAG 

IFT_FLAG: switch for the generation of particle formation time. 

= 0: generate formation time. 

= 1: read in formation time from zpc.inp. 

IRE_FLAG: switch to read initial conditions from zpc.inp. 

= 0: read the initial conditions. 

= 1: do not read the initial conditions. 

IGEN_FLAG: switch to generate initial conditions. 

= 0: do not generate initial conditions. 

= 1: generate initial conditions. 

IBST_FLAG: choice of global frame. 

= 0: the global frame (the ordering frame) is the collider center of mass 
frame for Au on Au collisions. 

= 1: means the global frame is taken to be the target frame. 
COMMON /PARA5 / LCONFIG, IORD_SCH 

LCONFIG: choice of geometric configuration and space cell division 
optimization. 

= 1: the system is undergoing 3-d expansion. 

= 2: the system is undergoing 1-d expansion with space cell division. 
= 3: the system is undergoing 1-d expansion without space cell division. 
= 4: the system is confined in a box with space cell division. 
= 5: the system is confined in a box without space cell division. 
IORD_SCH: choice of collision scheme. IORD_SCH = 10 * II + 12 
II: choice of collision prescription. 

= 0: the collision frame is the two parton center of mass frame. The 
scattering space point is the center point of the 2 parton positions at 
closest approach in the collision frame. 

= 1: the collision frame is the two parton center of mass frame. The 
scattering point is the position of the parton at the closest approach 
in the collision frame. This gives two generally different collision times 
for the two colliding partons. 
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= 2: the collision frame is the same as the ordering global frame. 
12: choice of ordering time. 

= 0: the earlier of the two collision times in the global frame is the 
ordering time for parton collisions. 

= 1: the average of the two collision times is taken as the ordering 
time. 

= 2: the later collision time is the ordering time. 

• COMMON /PARA6 / CENT_RAP 

CENT_RAP: the central rapidity of the rapidity plateau. 

• COMMON /PARA7/ LOSCAR 

LOSCAR: option to write OSCAR standard output. 

= 0: write OSCAR standard output. 

= 1: do not write OSCAR standard output. 

Six particle record common blocks are labeled prec!, in which ! goes from 
1 to 6. 



• COMMON /PREC1/ 

Purpose: stores the initial conditions for one event. 

• COMMON /PREC2/ 

Purpose: main particle record that has particle types, the positions 
at formation times, the corresponding formation times and the four 
momenta of the particles, plus the masses of the particles. 

• COMMON /PREC3/ 

Purpose: the same as /PREC2/ except the two colliding particles don't 
update their records until next collision. 

• COMMON /PREC4/ 

Purpose: stores velocities of particles. 

• COMMON /PREC5/ 

Purpose: stores formation time pseudo-rapidities, rapidities and proper 
times of particles. 
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• COMMON /PREC6/ 

Purpose: the same as /PREC5/ except the two colliding particles don't 
update their records until next collision. 

Six interaction list common blocks Llist! give us particle run time infor- 
mation and cell information. 

• COMMON /LLIST1/ ISCAT, JSCAT, NEXT(NMAXGL), LAST(NMAXGL), 
ICTYPE, ICSTA(NMAXGL) , NIC(NMAXGL), ICELSTA(NMAXGL) 
ISCAT: particle index. If the operation involves only one particle, then 

iscat is the particle index; if it involves 2 particles, iscat is the larger 
particle index. 

JSCAT: particle index. If the operation involves only one particle, then 

jscat is 0; if it involves 2 particles, jscat is the smaller particle index. 

NEXT(I): the next operation partner of particle i. 

LAST (I): the last operation partner of particle i. 

ICTYPE: the operation type. 

= 0: a collision between particles. 

= 1: the formation of a particle. 

= 2: both a collision between particles and the formation of a particle. 
= 3: a collision with wall. 

= 4: both a collision between particles and a wall collision. 
= 5: both a wall collision and a formation. 

= 6: a formation, collision between particles and a wall collision at the 
same time. 

ICSTA(I): the operation type for particle i. 

= 0: an ordinary collision. 

= 101: a collision with the wall with larger x. 

= 102: a collision with the wall with smaller x. 

= 103: a collision with the wall with larger y. 

= 104: a collision with the wall with smaller y. 

= 105: a collision with the wall with larger z. 

= 106: a collision with the wall with smaller z. 

= 111: a collision with another particle and the wall with larger x. 

= 112: a collision with another particle and the wall with smaller x. 

= 113: a collision with another particle and the wall with larger y. 

= 114: a collision with another particle and the wall with smaller y. 
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= 115: a collision with another particle and the wall with larger z. 
= 116: a collision with another particle and the wall with smaller z. 
NIC(I): the next particle index in the same cell as i. 
ICELSTA(I): the encoded information of the cell number particle i is 
in. If particle is in cell (il, i2, i3), then it equals il x 10000-M2 x 100-M3 
for particle inside the cube. When a particle is outside the 10 by 10 by 
10 box, its value is 111111. 

COMMON /LLIST2/ ICELL, ICEL(10, 10,10) 

ICELL: pointer to one particle out side the 10 by 10 by 10 cell volume. 
ICEL(I, J, K): one particle in the cell (i, j, k). 

COMMON /LLIST3/ SIZE1, SIZE2, SIZE3, VI, V2, V3, SIZE 

SIZE1: cell size in the x direction. It must be larger than the inverse 

square root of cutoff2. 

SIZE2: cell size in the y direction. 

SIZE3: cell size in the z direction. 

VI: cell expanding velocity in x direction. 

V2: cell expanding velocity in y direction. 

V3: cell expanding velocity in z direction. 

SIZE: the time that cells begin to expand. 

COMMON /LLIST4/ IFMPT, ICHKPT, INDX(NMAXGL) 
IFMPT: the particle index of the parton that is to be formed next. 
ICHKPT: the last formed particle index. 

INDX(I): the index of particle i of /INLREC/ in an increasing forma- 
tion time order. 

COMMON /LLIST5 / CT(NMAXGL), OT(NMAXGL), TLARGE 
CT(I): the collision time for particle i. 
OT(I): the ordering time for particle i. 
TLARGE: a large number for time cutoff. 

COMMON /LLIST6/ LOPERATION, LCOLLISION, T 

N_OPERATION: the current number of operations. 

LCOLLISION: the current number of collision between particles that 

have been performed. 

T: the current operation time. 
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4.4 Other common blocks 

Two common blocks /AUX_REC1/ and /AUX_REC2/ contain auxiliary vari- 
ables for 1-d expansion and particles in a box without space cell division. 
Two common blocks contain information about the random numbers. 

• COMMON /RANDOM1/ NUMBER 

NUMBER: the number of random numbers that have been used. 

• COMMON /RANDOM2/ IFF 

IFF: choice of attractive or repulsive force, it alternates between 1 and 
-1. 

Currently, there are four common blocks ana! for analysis. 

• COMMON /ANA1/ TS(12) 

TS: stores the time points that we want to sample data for analysis. 

• COMMON /ANA2/ 

Purpose: to record dE T /dy, dE T /dN, and dN/dy time evolution. 

• COMMON /ANA3/ EM(4, 4, 12) 

EM: energy momentum tensor for 12 different time points. 

• COMMON /ANA4/ FDETDY(24), FDNDY(24), FDNDPT(12) 
FDETDY: final dE T /dy distribution. 

FDNDY: final dN/dy distribution. 
FDNDPT: final dN/dp T distribution. 

4.5 Input files 

• zpc.ini 

Purpose: provides the following parameters for ZPC: 
'mass of particles (GeV)' 
'screening mass (fm(-l))' 
'alpha' 

strong interaction coupling constant; 
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'rts_cut2 (GeV2)' 

square root s cutoff bellow which no collision will by taken into account; 

'n_sta_evt' 
starting event number; 

'n_event' 
number of events; 

'n_run' 
number of runs; 

'gene ft flag (O,yes;other,no)'; 

'read ini flag (0, read; other, no)'; 

'gene ini flag (0,no;other,yes)'; 

'boost ini condition (0,no)'; 

'Lconfig (1,3-d exp;2,3 1-d exp;4,5 box)' 

'collision ordering scheme' gives the parameter IORD_SCH; 

'optional OSCAR output: (0,no;l,yes)' 

'vl, v2, v3' 
vl, v2, v3 are the cell expanding velocities; 

'sizel, size2, size3 (fm)' 
sizel, size2, and size3 give the cell sizes in the x, y, z directions; 

'beginning iff' 

beginning iff gives the beginning random choice of attractive or repul- 
sive collision; 

'iseed' 

iseed is the random number seed; 
'i_ran_used' 

number of random numbers used. This makes it possible to continue a 
series of events. 

• zpc.inp 

Purpose: provides initial conditions for ZPC Each particle has one 
record in zpc.inp. Each record has 11 entries: the event number, the 
particle type (integer), the x, y, and z positions of the particle at the 
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time specified, the specified time, and the four momentum of the par- 
ticle, the mass of particle (double precision). The specified time can 
be one time same for all the particles, in this case, the ift_flag in the 
zpc.ini file should be set to 0. The specified time can also be the given 
formation time of the particle. In this case, the ift_flag should be set 
to 1. 



4.6 Output files 

• zpc.res 

zpc.res is the major output file. It has the event number, run number, 
number of operations, number of collisions between particles, freezeout 
time, ending random number, and ending randomizer iff. The user can 
dump other analysis information into zpc.res. 

• zpc.jun 

zpc.jun saves execution time and error message. 

• zpc.oscar 



ZPC 1.0.1 writes OSCAR standard output |T2| into zpc.oscar. 



anal directory 

Currently files in anal are for the Lorentz invariance analysis, fdetdy.dat 
is the final detdy distribution; fdndy.dat is the final dndy distribution; 
fdndpt.dat is the final dndpt distribution. 

ana2 directory 

Currently files in ana2 are for the time evolution of detdy dndy and 
detdy dn analysis, detdy*. dat, dndy*.dat, and detdydn*.dat are for the 
time evolution of detdy, dndy, detdydn correspondingly. * is a number 
that specifies the rapidity range of the central slice that the program 
takes data from. 

ana3 directory 

ana3 has information for 1-d expansion semi-analytic tests, coll.res has 
the collision proper time and Mandelstam s for central collisions. et*.res 
files give lab time development of dEx/dy at fixed proper times. et*.res 
files are useful in choosing a final freezeout time for 1-d expansion such 
that the dE T / dy\ y=0 (r) reaches asymptotic value. 
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• ana4 directory 

ana4 contains energy-momentum tensor at different times for particles 
in a box case. 

7 Other files 

• README 

README gives a brief description of the files of ZPC. 

• COPYING 

COPYING contains GNU General Public License. 

• zpc.doc 

zpc.doc has detailed discussion of the program and instruction on run- 
ning the program. 

• zpc.f 

zpc.f is the main source file. 

• zpc.go 

zpc.go is the script file that controls the running of the program. On 
some machines, there are special rules for batch jobs, zpc.go must be 
modified accordingly. 

• zpc.mk 

zpc.mk is the makefile that generates the executable file. 

Instructions for users 

• Set parameters in zpc.ini; 

• Make zpc.inp if the initial conditions are to be read in from it; 

• Run the program by running zpc.go. It will use zpc.mk to make the 
executable file, run and time the executable file. 

Here is an example input file and corresponding results. 
Input file zpc.ini: 
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'mass of particles (GeV) ' OdO 

'screening mass (fm~(-l)) ' 3d0 

'alpha ' 0.47140452 

'rts_cut2 (GeV"2) ' O.OldO 

'n_sta_evt ' 1 

'n_event ' 10 

' n_run ' 1 

'gene ft flag (0,yes; other, no) ' 1 

'read ini flag (0,read;other,no) ' 1 

'gene ini flag (0 , no ; other , yes) ' 1 

'boost ini flag (0,no;other,yes) ' 

'i_config (1,3-d exp;2,3 1-d exp;4,5 box) ' 1 

'collision ordering scheme ' 11 

'optional OSCAR output: (0,no;l,yes) ' 

'vl, v2, v3 ' 0.2d0 0.2d0 0.2d0 

'sizel, size2, size3 (fm) ' 1.5 1.5 0.7 

'beginning iff ' -1 

'iseed ' 7 

'i_ran_used ' 1 



Output file zpc.res: 

Event 1 , run 1 , 

number of operations = 36224, 

number of collision between particles = 6914, 

freezeout time= 70119.524722901, 

ending at the 45159th random number, 

ending collision iff= -1 

Event 2 , run 1 , 

number of operations = 35886, 

number of collision between particles = 6856, 

freezeout time= 447736.85516922, 

ending at the 90145th random number, 

ending collision iff= -1 

Event 3 , run 1 , 

number of operations = 36118, 
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Figure 1: 10 events, [i — 3 fm 1 

number of collision between particles = 6954, 
freezeout time= 99267.891489987, 
ending at the 135343th random number, 
ending collision iff= -1 



The data files in directory anal give the final -r^ 1 , ^ and f^-. They are 

J ° ay 1 ay apx J 

shown in Fig. 1, Fig. 2 and Fig. 3. 

6 Modification and further development 

The current version 1.0.1 has only elastic collisions. It can be easily modified 
to study more complicated physical process, for example, to study effects of 
elastic collisions on chemical equilibrium. To put in more particle species: 

• Add in new parameters in zpc.ini and change subroutine ini_zpc corre- 
spondingly. 
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Figure 2: 10 events, \x = 3 fm 1 
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Figure 3: 10 events, ji = 3 fm 1 
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• Add in zpc.inp more entries corresponding to new degrees of freedom 
for particles and change read_ini correspondingly. 

• In subroutine scat, new subroutine is necessary to update the new 
degrees of freedom for the colliding particles, subroutine get_that needs 
to be changed correspondingly to get the scattering angle through the 
new scattering differential cross section. 

• In subroutine iscoll*, new subroutine get_cutoff2 is necessary to get - 
for different processes. 

Further analysis can be done by putting the analysis in the subroutine 
zpc_anal and zpc_ana2, and adding corresponding output in zpc_out. 
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